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A theoretical analysis is presented of a nematic liquid crystal confined between substrates pat- 
terned with squares that promote vertical and planar alignment. Two approaches are used to eluci- 
date the behavior across a wide range of length scales: Monte Carlo simulation of hard particles and 
Frank-Oseen continuum theory. Both approaches predict bistable degenerate azimuthal alignment 
in the bulk along the edges of the squares; the continuum calculation additionally reveals the possi- 
bility of an anchoring transition to diagonal alignment if the polar anchoring energy associated with 
the pattern is sufficiently weak. Unlike the striped systems previously analyzed, the Monte Carlo 
simulations suggest that there is no "bridging" transition for sufficiently thin cells. The extent to 
which these geometrically patterned systems resemble topographically patterned substrates, such 
as square wells, is also discussed. 



I. INTRODUCTION 

The imposition of a liquid crystal's (LC's) bulk direc- 
tor orientation through that LC's interaction with a con- 
fining substrate is termed anchoring [T]. In the absence 
of defects and applied fields, substrate anchoring is the 
main determinant of the director profile in a sandwich 
geometry LC cell; the director profile in such a cell is set 
through of the orientational elastic energy, subject to 
each wall's polar and azimuthal anchoring constraints. 
At a continuum level, this elastic energy is most com- 
monly expressed through square director gradient terms 
corresponding to the independent splay, twist and bend 
modes of orientational deformation, weighted by the elas- 
tic constants K\, K2 and K$, respectively. When consid- 
ered at a finer length-scale, conversely, both the anchor- 
ing strengths and the bulk elastic constants are emer- 
gent from the microscopic interplay of the orientational 
and positional degrees of freedom of the liquid crystalline 
molecules and the confining surfaces. 

Traditional routes to establishing desired anchoring 0, 
and, thus, director profiles of use for display devices, 
include substrate rubbing and various photo-alignment 
approaches (light-induced cis-trans isomerisation, pho- 
todegradation ...). Whilst significant empirical knowl- 
edge has been developed in relation to each of these 
approaches, the molecular mechanisms by which they 
achieve desired are poorly understood. Introducing 
inhomogeneity into substrate conditions has, for some 
time, been as an alternative route to both control- 
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ling conventional anchoring and, increasingly, introduc- 
ing new phenomena. A range of such substrates have 
been developed and examined. These cover patterning 
length-scales ranging from 10 m upwards, couple to the 
LC either sterically, chemically or dielectrically (or by a 
combination of same) and have been achieved as both 1- 
dimensional (stripes, ridges, ...) and 2-dimensional (cir- 
cles, squares, triangles, posts, ... ) patternings. 

One of the important phenomena that can be achieved 
through substrate patterning is bistability, that is of 
two distinct anchoring arrangements with (in the absence 
of an applied field) mutually inaccessible free energy min- 
ima. Pattern-0 has now been established for the blazed 
grating structure of the Zenithally Bistable Device [5J [3] , 
the two-dimensional array of the Post Aligned Bistable 
Nematic |4] device and, more recently, a steric square- 
well arrangement [Sl[6|. In each of these, the bistability 
pertains between one state with a continuous director 
arrangement and a second containing orientational de- 
fects that are pinned in some way by the substrate in- 
homogeneity. This suggests that the key length-scale for 
achieving bistability here is the size and periodicity of the 
patterning, a conjecture which is supported by the suc- 
cess of mesoscopic modelling approaches in both access- 
ing the bistable states and relating, semi-quantitatively, 
switching fields to geometrical parameters pj |5] ■ 

In addition to these sterically patterned systems, 
chemical patterning has now also been developed as an 
approach for imposing substrate inhomogeneity on LC 
systems. The notion of imposing combinations of az- 
imuthal and polar anchorings on LCs via chemically 
nanopatterned substrates was the subject of early ex- 
perimental work [5HTT]. Subsequently, Lee and Clark 
performed a more systematic study of the alignment 
properties of nematic LCs on surfaces comprising both 
homeotropic and planar alignment areas |12j . For stripe 
patterns, they found that the polar orientation depends 
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on the relative areas of the homeotropic and planar re- 
gions, but that the azimuthal anchoring always runs 
along the direction of the stripes. In 2005, Scharf et 
al. [131 Q3] undertook further investigations of systems 
with competing alignment regions. Further innovations 
by the groups of Abbott [15] and Evans [16] on the 
development of patterns of combinations of alkanethiols 
deposited as self assembled monolayers (SAMs) on gold. 
Using microcontact printing, these systems proved ca- 
pable of achieving highly reproducible surface features 
with periodicities of 10s of microns. Square, circular and 
striped patterns written on these length-scales were, thus, 
observed using optical microscopy in crossed 2 set-ups. 
An alternative approach, employing selective ultra-violet 
irradiation of SAMs, achieved LC-aligning stripe patterns 
on the sub-micron scale |17j . 

In Ref [18] , two of the current authors contributed the 
simulation aspects of a joint experimental / simulation 
study of LC alignment at a single patterned substrate. 
In this, it was shown that a range of patterned SAMs can 
be used to control LC alignment states and domains. For 
stripe patterns, the LC was found to align parallel to the 
stripe boundaries for both nanoscale simulation features 
and micron-scale experimental systems. Indeed, despite 
the significantly different length-scales involved, the qual- 
itative seen in simulations of generic molecular models 
confined using a striped substrate proved entirely con- 
sistent with the experimental observations. Specifically, 
on undergoing isotropic to nematic ordering, all systems 
proved to be dominated by the homeotropic-aligning sub- 
strate regions at the ordering transition, the influence of 
the planar-aligning regions only becoming apparent well 
into the nematic phase. 

In Ref |19) , we extended our molecular simulation work 
to consider the of a thin nematic film confined be- 
tween two identical nano-patterned substrates. Using 
patterns involving alternating stripes of homeotropic-0 
and homogeneous-0 substrate, we showed that the po- 
lar anchoring angle can be varied continuously from pla- 
nar to homeotropic by appropriate tuning of the rela- 
tive stripe widths and the film thickness. For thin films 
with equal stripe widths, we also observed orientational 
bridging, the surface patterning being written in domains 
which traversed the nematic film. This dual-bridging- 
domain arrangement broke down with increase in film 
thickness, however, being replaced by a single tilted mon- 
odomain. Strong azimuthal anchoring in the plane of the 
stripe boundaries was observed for all systems. 

Stripe-geometry systems have also been analyzed by 
the third of the current authors using continuum the- 
ory [20 22 . This larger length-scale work has shown that 
the basis for azimuthal alignment by striped substrates is 
associated with differences in the Frank elastic constants. 



Azimuthal anchoring parallel to the stripes corresponds 
to the LC adopting a configuration comprising twist, 
splay and bend deformations; in the other limiting case, 
bulk alignment perpendicular to the stripes, only splay 
and bend deformations are required. Experimentally, Ki 
is significantly lower than K\ and for most nematics, 
so that parallel anchoring is stable. Monte Carlo esti- 
mates of the elastic constants for calamitic particle-based 
LC simulation models yield similar elastic constant ra- 
tios [23] , so this phenomenological agreement between the 
predictions of particle-based and continuum approaches 
is to be expected. 

In this paper, we extend our respective works on stripe- 
patterned systems by studying the effect of substrates 
with square patternings on a confined LC film. Experi- 
mental studies of such systems include the checkerboard 
patternings achieved by Bramble [18] and, more recently, 
Yi [24] and the bistable square-well systems mentioned 
above [3[B]. In respect to the latter, we note that both 
Q-tensor [55] and Landau-De Gennes [5f5] modelling ap- 
proaches have been used to examine the stable configu- 
rations for such systems. From this, diagonally-anchored 
and edge-anchored states have been identified the former 
comprising surface region defects. 

Here we use both molecular- and continuum-level mod- 
elling approaches to investigate the of LC films con- 
fined between square-patterned substrates. In section 
II we present our molecular-level model system and de- 
scribe the simulation methodology employed. Section III 
contains the corresponding simulation results. Following 
this, in Section IV we present a continuum- level analy- 
sis of anchoring control in systems with square-patterned 
substrates. Finally, in Section V, we compare and com- 
bine the findings from these investigations to draw more 
general conclusions. 

II. MOLECULAR MODEL AND SIMULATION 
DETAILS 

We have performed a series of Monte Carlo (MC) sim- 
ulations of rod-shaped particles confined in slab geome- 
try between two planar walls. Inter-particle interactions 
have been modelled through the Hard Gaussian Overlap 
(HGO) potential [57]. Here, the dependence of the inter- 
action potential v HGO on U; and Uj, the orientations of 
particles i and j, and fy, the inter-particle unit vector is 

{0 if nj > a(rij,Ui,Uj) 
oo if nj < a(iij,Ui,Uj) (1) 

where <r(ry, Uj, u,-), the contact distance, is given by 
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<r(f^,Ui,Uj-) = cr 
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1 + x(ui.Uj) 
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The parameter % is set by the particle length to breadth 
ratio k= a en d /a S ide via 
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Particle-substrate interactions have been modelled us- 
ing the hard needle-wall potential (HNW) In this, 
the particles do not interact directly with the surfaces. 
Rather the surface interaction is achieved by considering 
a hard axial needle of length aok s placed at the center of 
each particle (see Figure [I]). This gives an interaction 



,HNW 



if \zi - z \ > a w (ui) 

oo if \zi - z \ < a w (ui) (4) 



where zq represents the location of a substrate and 
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distance L z = 4Ker , periodic boundary conditions being 
imposed in the x- and y-directions. 

On each substrate, k s was set to a homeotropic- 
aligning value (fc s = 0) for two quadrants of its area and 
a planar value (k s = 3) for the remainder. Except where 
otherwise stated, sharp boundaries have been imposed 
between the different alignment regions. The patterns 
on the top and bottom surfaces have been kept in perfect 
registry with one another, as shown in the schematic in 
Fig. [2j Each system was at low density and compressed, 
in small increments, by decreasing the box dimensions L x 
and L y whilst maintaining the condition L x /L y = 1. At 
each density, a run length of 1 million MC sweeps (where 
one sweep represents one attempted move per particle) 
was performed, averages and profiles being accumulated 
for the final 500 000 sweeps. 




FIG. 2: (Color online) Schematic representation of rectan- 
gle patterned systems with alternating homeotropic-inducing 
(red/dark) and planar-inducing (green/light) substrate re- 
gions. The Euler angle is from the x-axis. 



FIG. 1: (Color online) Schematic representation of the geom- 
etry used for the hard needle-wall (HNW) particle-substrate 
interaction [28|. 

arccos(u; i2 ) is the angle between the substrate normal 
and the particle's orientation vector. For small k s , the 
homeotropic arrangement has been shown to be stable, 
whereas planar anchoring is for long k s [28!. Fur- 
thermore, despite its simplicity, the HNW potential has 
been found to exhibit qualitatively identical to that ob- 
tained using more complex particle-substrate potentials 
[29) . Here, by imposing variation in k s across the two 
boundary walls, we investigate the effects of molecular- 
scale substrate patterning on LC anchoring. The results 
presented in Section [Hl| were obtained for systems of 864 
k = 3 HGO particles confined between two square pat- 
terned substrates. The substrates were separated by a 



Analysis was performed by dividing stored system con- 
figurations into 100 equidistant constant-z slices and cal- 
culating averages of relevant observables in each slice. 
This yielded profiles of quantities such as number den- 
sity, p*{z), from which structural changes could be as- 
sessed. Oricntational order profiles were also measured, 
particularly 




which measures variation across the confined films of ori- 
entational order measured with respect to the substrate 
normal. Here N(z) is the instantaneous occupancy of 
the relevant slice. We have also further subdivided the 
system to assess lateral inhomogeneities induced by the 
patterning. 
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III. MONTE CARLO SIMULATION RESULTS 

In the initial system considered here, sharp transitions 
were imposed between the k s = and k s = 3 regions. 
The outcomes of these squared system simulation are 
by the snapshots shown in Fig.[3j Several remarks emerge 
from these. The substrate patterning is readily apparent 



this system more quantitatively by calculating two sets 
of profiles of key observables for this system; for anal- 
ysis purposes, each simulated system has been split in 
two according to the imposed substrate pattern. In this, 
individual particles have been allocated to homeotropic- 
confined or homogeneous-confined regions according to 
their x- and y-coordinates. 




(d)p*=0.38 



(e)p*=0.40 



FIG. 3: (Color online) Snapshots of the square patterned sys- 
tem with sharp transitions between k s =0 and k a —3 regions 
for a series of different reduced densities. Particles are coded 
for orientation: x red; y green; z blue. 

from all of these, with ordered layers of homeotropic- 
aligned and planar-aligned particles residing in the ap- 
propriate regions. Sharp delineation between these re- 
gions can be seen for all densities. At p*=0.30 and 0.34, 
the particles at the of the film appear to be relatively 
disordered, whereas aligned monodomains can be seen at 
the three higher densities. Animations of these simula- 
tions show that in the planar aligning substrate regions, 
the molecules regularly flipped en masse between the x- 
and y- orientations. This tendency is apparent by com- 
paring Figs. 3(c)| |3(d")| and 3(e) the orientations on the 
planar part of the substrate vary from image to image. 
At high density (p*=0.40), the planar parts of the sur- 
face patterns appear restricted to the immediate vicinity 
of the substrate regions. Unlike the equivalent stripe pat- 
terned svstem[T§]. for which the patterning fully bridges 
a film of this width, the planar parts of the pattern ap- 
pear unable to penetrate the LC film. With this very 
short length-scale pattern, then, the spontaneous flipping 
of the planar substrate domains between the two degener- 
ate arrangements appears to make it impossible to fix or 
control the azimuthal anchoring orientation. At p*=0.37, 
the system appears to exhibit homeotropic anchoring 
(Fig. |3(c)| ). On increasing the density (p* =0.38), how- 
ever, this initial homeotropic anchoring adopts a clear tilt 
(Fig. |3(d)[ ) through which the planar alignment regions 
on the two substrates become coupled. On further com- 
pression to p*=0.40, the bulk director parti ally r egains 
its alignment normal to the substrates (Fig. |3(e) ). 
In the light of these observations, we have the of 
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(a)Density profile: homeotropic-confined region 
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(b)Density profile: planar-confined region 

FIG. 4: (Color online)Density profiles for the square 
terned system at different reduced 1 p*. 



pat- 



The density profiles depicted in Fig. 4(a) show the ad- 
sorption characteristics for the portion of the film con- 
fined between the homeotropic surface regions. These 
indicate that increasing the density leads to formation of 
surface layers with a period icity of ~2<7o (i.e. 2/3 of the 
particle length). Fig. |4(b)| shows the corresponding of 
the regions of the film confined between the planar- con- 
fining surfaces. Here, a shorter wavelength density oscil- 
lation is apparent close to the substrates. Despite these 
differences close to the substrates, both profiles adopt 
very similar in the central part of the film: essentially 
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featureless at low (isotropic) densities and oscillatory at 
high (nematic) densities. These oscillations are consis- 
tent with the formation of a homeotropic (or near 0) bulk 
monodomain. Such formation is only seen for much 
thicker films when stripe patterning is imposed |19j . 

A more complete understanding of the orientational 
aspects of the substrate-induced ordering in this system 
can be obtained from the Q zz diagonal component of the 
order tensor. For perfect homeotropic anchoring, Q zz {z) 
should tend to 1 and for perfect planar anchoring, Q zz (z) 
should tend to -0.5. Fig. 5(a) shows the Q zz profiles mea- 



sured in the homeotropic-confined regions. As the den- 
sity is increased, initially the bulk-region Q zz value in- 
creases as well, showing the development of homeotropic 
anchoring in the bulk. At a density of 0.37, the bulk 
Q zz value reaches 0.60-0.65. On further increasing the 
density to 0.38, however, the Q zz value decreases to just 
below 0.5. Then, as the density reaches 0.4, the Q zz in- 
creases again to Q 2 z=0.60-0.65. This non-monotonic 
confirms, in a statistically significant fashion, the tilt 
apparent in the corresponding snapshots. 
It is also informative to compare these observations 
with equivalent profiles obtained for HGO films con- 
fined between unpatterned homeotropic- and planar- 
aligning substrates. To this end, Fig. [6] shows that 
for the equivalent unpatterned homeotropic-aligning sys- 
tem, increasing the density causes the central Q zz value 
to increase monotonically(red line-green line-dark blue 
line/gray line-light gray line-dark gray line). The fact 
that Q zz shows a decreases at p* ~0.37 in the patterned 
system is, then, associated with a tilt of the bulk director 
caused by the planar pattern regions on the surface. 

An equivalent comparison performed for the planar- 
aligning region (Fig. [8]) shows a marked difference between 
the patterned and unpatterned systems. Indeed, despite 
its intrinsic anchoring 0, the Q zz (z) of the planar- 
aligned region of the patterned system is actually far 
closer to that of the unpatterned homeotropic-confined 
system. Only very close to the substrates is the planar 
nature of the imposed substrate pattern apparent. 

In order to assess the azimuthal anchoring in this 
system, we also constructed a time-averaged histogram 
of the molecular angles observed during the p*=0.37 
simulation. Specifically, this histogram (Fig . [7]) was 
generated from 500 stored configuration files and based 
on the orientations of particles within Ictq of the pla- 
nar substrate regions. The histogram is strongly peaked 
at angles corresponding to the boundaries of the square 
pattern, i.e. the molecules at the planar substrates 
are strongly disposed to adopting azimuthal angles cj> of 
^ 0° < <f> < 5° and ~ 90° < <j) < 90°. This is consis- 
tent with our previous observation that the molecules on 
this region appeared to regularly flip between the x- and 
y-directions. 

Before closing this section, we return to the observation 
that the substrate patterning applied here failed to pene- 
trate the LC film beyond the first adsorbed monolayer. In 
fact, we can report that this was a general characteristic 
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(a)Q zz profile for the homeotropic-confined region 



O 




(b)Q zz profile for the planar-confined region 

FIG. 5: (Color online)Q zz profiles for the square patterned 
system at different reduced 1 p* . 



observed for a range of different two-dimensional pattern- 
ings; simulations performed with circle, oval and rectan- 
gle patterns, respectively, all led to the development of 
central monodomain configurations. To explain why the 
orientational bridging observed in stripe-patterned sys- 
tems is lost on moving to two-dimensional patterns, we 
consider the hypothetical square bridging domain shown 
schematically in Fig. [9j Here a planar-aligned domain is 
bounded at each face by homeotropic material. From the 
schematic, though, it is clear that two distinct pairs of do- 
main boundaries are required for this: one pair involving 
T-like orientational changes and another involving X-like 
configurations. Such a scenario is clearly unstable: the 
surface tensions associated with these different pairs of 
interfaces could not, in general, be equal. As a conse- 
quence, the hypothetical square bridging domain consid- 
ered here could never be a stable arrangement. Indeed, 
similar instabilities hold for all orientational bridges pro- 
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FIG. 6: (Color online)Comparison of the Q zz profiles for 
unpatterned (lines) and square-patterned (symbols) systems: 
homeotropic regions. 



FIG. 8: (Color online) Comparison of the Q zz profiles for 
unpatterned (lines) and square-patterned (symbols) systems: 
planar regions. 
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FIG. 7: (Color online)Histogram showing the azimuthal angle 
distribution of surface particles in the P region. 



jected from two-dimensional patternings. 

To conclude, these simulations indicate that LC be- 
tween square-patterned substrates have a tendency to 
form monodomains. These monodomains are different 
from those developed between unpatterned substrates, 
though, since a) they can exhibit a non-monotonic 
density-dependent tilt and b) the azimuthal anchoring 
shows a strong coupling parallel to the square edges but 
is degenerate between the different edge orientations. 



IV. CONTINUUM MODEL 

To further understand the aligning effects observed in 
the simulations presented in the previous section, we now 
consider the of systems with the same geometry of pat- 



T configuration 




X configuration 



FIG. 9: (Color online)Schematic x-y plane slice through a 
hypothetical planar aligned bridging domain in a square pat- 
terned system. T-configuration : splay/bend distortion. X- 
configuration : Twist 



terning but applied at a much larger length-scale. Specifi- 
cally, we analyze the effect of square-patterned substrates 
on LC films in the continuum limit. In this approach, the 
local orientation of the nematic is characterized by a unit 
vector field known as the director and here by 



n(a ? ) = (cos 9 sin cf>, cos 9 cos (f>, sin 9) 



(7) 
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where the coordinates are chosen as depicted in fig. 2. 
The actual configuration adopted by the nematic is that 



which minimizes the Frank free energy 



F=^Jd 3 xKi(V- hf + K 2 [n • (V x n)] 2 + K 3 |n x (V x n)| 2 + J AS g{n, n ). 



(8) 



r 



Here, the first integral is to be performed over the volume 
of the nematic layer and the second over the surfaces in 
contact with the 0. The interaction of the nematic with 
the surface is characterized by an anchoring potential 
g(ri, no) that measures the energy cost of moving the di- 
rector away from an easy axis no; for a patterned surface 
this varies as a function of position. The configuration of 
the LC is found by solving the Euler-Lagrange equations 
for 9 and <f>; these are generally nonlinear and difficult to 
solve analytically in more than one spatial dimension. 

A common simplification, known as the one-constant 
approximation, is to set K\ = K 2 = K 3 . If this is done, 
the Euler-Lagrange equation for 9 reduces to Laplace's 
equation. However, such an approximation is unsuitable 
for analyzing situations with patterned surfaces because 
the aligning effect on the LC is due to differences between 
the elastic constants [30]. A "two- const ant" approxima- 
tion, where K\ — K% ^ K% has been previously used 



to understand the situation of a nematic film in contact 
with a surface patterned with alternating homeotropic 
and planar stripes, a two-dimensional system [20], [22] and 
here we extend the analysis to three dimensions. 

In order to proceed, a further simplifying assumption 
is made: that the director is confined everywhere to a 
single plane, i.e. that cj> is spatially uniform. This sim- 
plification is motivated (and justified) by the observation 
that the molecular distribution of azimuthal angles (Fig 
. [7| in our MC simulations implies a monodomain ar- 
rangement for all nematic films. Whilst is taken to be 
constant, the zenithal, or tilt, angle 9 remains free to vary 
in response to the substrate pattern. Keepingp/iifixed is 
further motivated by the observation that wherever the 
director is nearly homeotropic, variations in <p contribute 
negligibly to the free energy. The free energy density in 
this situation is 



TCOS 



+(1 - t) sin(20) 



89 

dx 

89 89 
dx 8y 



+ (rsii 



cos 



(9) 



and the corresponding EulerLagrange equation for 9 is 
linear 



tcos ffl + sm 
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This can be converted to Laplace's equation in new co- 
ordinates (£,f7,C) by the following linear transformation 



where 



x 

P T QP ( y 

z 



(11) 
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(12) 
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and 
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(13) 



The geometric interpretation of the transformation is a 
combination of a rotation and shear. To solve the Euler- 



Lagrange equation for 9(x, y, z), eq. 10 we try a solution 
of the form 
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9(x,y,z)=9 Q + J2 E x( A 
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+ B nm e 1, ' lmZ ) exp [i2n (nx + my) /A] 



r 



(14) 



where A = L x /L z such that XL Z is the period the pattern- if the parameters v nm arc chosen thus 
ing in both x and y directions. The equation is satisfied 



v nm = 7iV2(r + 1) (m 2 + n 2 ) - 2(r - 1) [2mn sm(2<j>) + cos(20) (to 2 - n 2 )] 



r 



(15) 



The constant 9q is, from the mean-value theorem, 
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(16) 



The coefficients A nm and B nm are determined by the 
boundary conditions. For weak anchoring, these are the 
torque-balance equation 



g df(9,V9) dg(9,9 e ) =Q 



dV9 
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(17) 



evaluated at each surface, where s is the outward sur- 
face normal. To facilitate separation of the in[l6j the 
harmonic anchoring potential 



g H (0,9 e ) 



2 



(18) 



is chosen, yielding the Robin boundary condition 

09 „ n 



Oz 



(19) 



where the — ve sign corresponds to z — z$ = —L z /2, 
the +ve sign to z — zq — +L z /2 and the dimensionless 
parameter associated with polar anchoring Lg is 



W e d' 



(20) 



Inserting the solution ( 14 ) into the boundary condition 



(191 at each surface yields the coupled system of equa- 



tions 



1 + L e v n m 1 — L s v nm 
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where c nm and d nm are the Fourier coefficients of the 
easy axis profile 9q(x, y) at the z — and z = 1 surfaces 
respectively. These are simply 



— dn.vn. 
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■Knra 



n, to odd 
otherwise 
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and solution of (21 ) yields 
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The complete director profile for given values of 4>, r, Lg 
and A is then fully specified by the series solution ( 14 ) 



and the parameters ( 15 1 and ( 23 ) that have now been 
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determined. 



The free energy associated with the solution ( 14 1 may 
be evaluated by substituting it into the free energy ^ 



and performing necessary integrations. The bulk energy 
is 



7T 

F b = V27; — [{ A nm e ~ Vnm + B nm e+Vnm ) smh{v nm ) + 2A nm B nm U nrn \ 



x {(1 + r) (m 2 + n 2 ) + (1 - r) [cos(20) (m 2 - n 2 ) + 2mnsin(2(/>)] } 



si nil ( i^nm ) 2 A nrn B nrn v n 



(24) 



r 



The surface energy (for each surface) is 

7t X /16-I- — y (A nrn -\-L> nm ')(A n . m -\-I3 nm 2c nm ). 

Lg — 

nm 

These expressions for the free energy have been eval- 
uated numerically as a function of r/> for different values 
of Lg. A value of r = K 2 /Ki = 1/2 was used that is 
approximately valid for many common nematics includ- 
ing 5CB. The period of the pattern was initially cho- 
sen to be the same as the cell thickness, i.e. A = 1. 
The plots displayed in fig. [TO] reveal an anchoring transi- 
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FIG. 10: Free energy of the nematic as a function of cf>, plotted 
for various values of Lg. 



tion: as Lg —¥ 0, representing rigid polar anchoring, the 
squares promote azimuthal alignment parallel to their 
sides, and there are two degenerate solutions at <j> — 
or <j) — 7r/2, i.e. the same as was observed in our 
molecular simulations. If Lg is increased, however, align- 
ment along the diagonals i.e. <fi = ±7r/4 becomes the 
energetically preferred solution. The critical Lg at which 
the diagonal and aligned solutions become degenerate is 
roughly Lg ~ 0.03. We can relate this result to experi- 
ment by taking K\ « 6pN (5CB), a typical cell thickness 



d — 10/im and a typical value for Wg = 10 -5 Jm~ 2 , which 
together give Lg = 0.006. This falls well into the regime 
of <j) — and tt/2 anchoring which, indeed, is precisely 
what was seen by Bramble [IB] . 

The second parameter of interest is A, the overall size 
of the squares relative to the cell thickness. Shown in fig- 
ure [TT] is the energy difference between the aligned = 
and diagonal 4> = 7r/4 solutions as a function of A, plot- 
ted for various values of Lg. In this plot, therefore, the 
diagonal solution is stable where lines lie below the ab- 
scissa, whereas the edge aligned solution is stable where 
the lines take positive values. We see, therefore, that 
the aligned solution is preferred both at small /lambda 
and as A — > 00. Reducing Lg has the effect of narrowing 
the window of A values for which the diagonal solution is 
preferred. 

For comparison, regarding the Monte Carlo simulation, 
A=1.12 at reduced density of 0.4 (Fig. 



3(e)) 
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FIG. 11: Free energy difference between the aligned = 
and diagonal <f> = 7r/4 solutions as a function of the period of 
the pattern A. 

Continuum predictions for Director tilt profiles as a 
function of z for the planar and homeotropic regions are 



displayed in fig. 12 and show that a nearly uniform con- 
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figuration is adopted at the cell center, similar to the 
findings from our particle-based simulations. As the an- 
choring parameter L g is reduced, the tilt conditions at the 
substrates relax and the central uniform region widens 
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FIG. 12: Calculated director angle profiles 6(z) in the center 
of the homeotropic and planar regions, for various values of 
Lg. 



CONCLUSIONS 



Alignment of a nematic between two substrates pat- 
terned with alternating homeotropic and planar squares 
has been studied using two theoretical approaches: A 
monte carlo simulation of rigid particles interacting 
through the hard gaussian overlap potential and a calcu- 
lation performed with nematic continuum theory. Both 
techniques show a regime where the nematic azimuthally 
aligns in the bulk with the edge of the squares. In the 
MC simulation, the average azimuthal orientation of the 
molecules is observed to flip between the two sides dur- 
ing the run; these states are energetically degenerate in 
the continuum approach. Furthermore, the continuum 
calculation reveals the existence of an anchoring transi- 
tion. If the polar anchoring is sufficiently weak and the 
period of the pattern is somewhat greater than the cell 
thickness, the nematic instead aligns along the diagonals. 
Unlike previously considered systems of square posts and 
wells 25, 26 our calculations seem to raise the possibil- 
ity of a re-entrant behavior as a function of the period of 
the pattern: with appropriate values of the polar anchor- 
ing energy, the diagonal state becomes unstable at both 
short and long wavelength patterning (fig. 11 ). It is likely 
that the location of the critical values of the period are 
only approximately correct because azimuthal variations 
of the director were not included in our minimization of 
the free energy. We expect that all of these transitions 



might be accessible with the Monte Carlo approach for 
larger simulation sizes; we are presently examining such 
systems. A transition between the two regimes ought to 
be experimentally observable by adjusting the ratio of 
the period of the pattern to the cell thickness. 



Our continuum analysis also revealed the 3 result that 
an exact linear form of the nematic Euler-Lagrange equa- 
tions exists even if there is 3D variation in the director 
and the nematic has elastic constants. The form of 
the resulting equation lends a geometric interpretation 
to elastic anisotropy as a transformation into skewed co- 
ordinates. This result should be of utility for further 
study of LC behavior in complicated geometries. 



The behavior observed in the present system is qual- 
itatively quite different to the case where a nematic is 
aligned by a striped surface. Previous Monte Carlo sim- 
ulations show that for a stripe system in a sufficiently 
film, there exists a regime where the nematic is divided 
into domains of vertical and planar alignment that bridge 
the corresponding regions of both substrates across the 
film. No similar behavior was observed for the square 
patterned system, where instead the planar part of the 
surface patterns is confined to the immediate vicinity of 
the substrate. This result may be interpreted by consid- 
ering the rotation of the molecular orientation between 
domains, which must occur through either a splay-like or 
twist-like distortion. Each of these types of wall would 
have a different energy density, described by a line ten- 
sion between domains with two different values. Since 
similar behavior appears from both monte carlo and con- 
tinuum theory, despite the fact that the latter entirely 
neglects variations in ordering, it appears that the line 
tension depends primarily on elastic distortion of the di- 
rector and does not significantly depend on the scalar 
order parameter. This justifies use of the continuum the- 
ory for such systems. 



In either alignment regime, whether along sides or di- 
agonals of the square pattern, the alignment is bistable 
due to the symmetry and hence of interest for electro- 
optic, display and sensing applications. In this re- 
spect these systems are quite similar to the post-aligned 
bistable display [35] or the arrays of square wells 
previously studied, which have the same symmetry, but 
there is an important difference: in the present case the 
bistable states have no disclinations present in the ne- 
matic configuration. Our results therefore reveal that 
defects are not necessary to stabilize bistable states. De- 
spite the apparent simplicity of the geometry, our results 
illustrate the rich phase diagrams exhibited by complex 
fluids in patterned geometries; further study of these sys- 
tems is presently being undertaken, to identify optimal 
switching strategies between the bistable states identified 
here. 
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